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ABSTRACT 

Variable ultraluminous X-ray sources (ULXs), which are considered to be black hole binaries (BHBs), 
are known to show state transitions similarly to Galactic BHBs. However, the relation between the 
ULX states and the Galactic BHB states is still unclear primarily due to less well-understood behaviors 
of ULXs in contrast to the Galactic BHBs. Here, we report a statistical X-ray spectral study of 34 
energy spectra from seven bright ULXs in the interacting galaxy systems M51 and NGC 4490/85, 
using archive data from multiple Chandra and XMM-Newton observations spanning for a few years. 
In order to compare with Galactic BHB states, we applied representative spectral models of BHBs; a 
power-law (PL), a multi-color disk black body (MCD), and a slim disk model to all the ULX spectra. 
We found a hint of a bimodal structure in the luminosity distribution of the samples, suggesting that 
ULXs have two states with typical luminosities of 3-6xl0 39 and 1.5-3xl0 39 ergs s _1 . Most spectra 
in the brighter state are explained by the MCD or the slim disk model, whereas those in the fainter 
state are explained by the PL model. In particular, the slim disk model successfully explains the 
observed spectral variations of NGC 4490/85 ULX-6 and ULX-8 by changes of the mass accretion rate 
to a black hole of an estimated mass of <40 Mq. From the best-fit model parameters of each state, 
we speculate that the brighter state in these two ULXs corresponds to the brightest state of Galactic 
BHBs, which is often called the "apparently standard state". The fainter state of the ULXs has a 
PL shaped spectrum, but the photon index range is much wider than that seen in any single state of 
Galactic BHBs. We thus speculate that it is a state unique to ULXs. Some sources show much fainter 
and steeper spectra than the faint state, which we identified as another state. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: individual (M51, 
NGC 4490, NGC 4485) — X-rays: binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are off-nuclear 
point-like sources detected in the X-ra y bandpass with 
lumin osi ties of >1 39 ~ 41 ergs s _1 ; see iLiu. fc Bregmanl 
(|2005h . iFabbianol (|2006t) for a review. Most vari- 
able ULXs are considered to be black hole binaries 
(BHBs) based on observational characteristics similar to 
Gal actic BHBs, such as sho rt- and long-term variations 
(e.g.lPtak. fc G riffiths 1999) , thermal disk em ission (e.g. 
iMakishima et al.l 120001 iFeng fc Kaaretll2010D . and state 
transitions fe.g. lKubota et alj|2001h . The large luminos- 
ity of ULXs suggests the black hole (BH) mass (M) to be 
much larger than that of Galactic BHs (~10 Mq) under 
the assumption that the radiation is spherically symmet- 
ric and the observed luminosity does not exceed the Ed- 
dington luminosity L^dd = 1.5x10 38 (M/Mq) ergs s -1 . 

Three competing interpretations have been proposed 
for ULXs. The first is that ULXs contain a so-called 
intermediate mass BH with a mass of 100-1000 Mq 
([Miller, fc Colbertl I2004D . In fact, a large mass of 
>700 Mq i s suggested for M 82 X-l based on it s larg e 
luminosity (|Matsumoto et al.ll200l1: iKaaret et~aT]|2001l ). 
The second is that ULXs are BHs with a mass compa- 
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rable to or slightly larger than tha t of Galactic BHs of 
<40 Mq (so-called s tellar mass BHs: lEbisawa et al.ll2003l: 
lOkaiima et al.l 2006) . The super-Eddington luminosity is 
inter preted as a consequence of the ULXs having a slim 
disk (lAbramowicz et allll988f) r ather than the standard 
disk (|Shakura. fc Sunvaevl ll973h In the state expressed 
by the slim disk model, which is a stable solution for very 
high mass accretion rates of >8x 10 18 (M/M Q ) g s _1 , 
the emission is moderately collimated toward the direc- 
tion normal to the di sk so that the lum inosity can be 
super-Eddington (e.g. lOhsuga et al1l2005t ). The third is 
the b eaming model with a high collimation (|King et al.l 
l2001h . although the scenario may have difficulties in ex- 
plaining extended photoio nized nebulae found around 
some ULXs. For example, Kaa ret et al.l (|2004h studied 
the Holmbergll ULX within an extended nebula. The 
nebula emits He II 4686 A line relatively isotropically, 
which is produced by the X-ray photoionization from the 
central ULX. Thus, the extended nebula is unlikely ex- 
cited by highly collimated radiation. 

Galactic BH Bs are known to show transitions among 
several states (|Esin et all I1997T ) . The most three well- 
established states are the low-hard state (LHS), the high- 
soft state (HSS), and the very high state (VHS). In 
general, the LHS is the faintest state (0.01-0.04 I/Edd 5 
iMaccaronel 120031 ) . which shows a pow er-law (PL) spec - 
trum with a photon index T = 1.5-1.9 (|Esin et al.lll997l) . 
The HSS is generally brighter than the LHS, and is 
a thermal state with a convex-shaped spectrum. The 
spectrum can be fitted by a multi-color disk black body 
(MCD) model (|Prmgldll98lh , which is an approximation 
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for the standard disk spectra. The VHS usually shows an 
X-ray luminosity comparable or higher than that of the 
HSS and a steep PL-like spectrum (typically f ~ 2.5, 
lEsin et ail Il997f) . The spectrum is considered to be a 
mixture of a weak disk component and a strong Comp- 
tonized component. 

In addition to these well-established three states, 
yet another state - - the app arently standard state 
(ASS) was proposed by iKubota. k Makishimal 

(2004). The ASS is a thermal state, which is 
brighter (Lx > 0.3 ^Edd) than the other states. At 
least two Galactic B HBs are reported to show this 



TABLE 1 
Observation Log. 



state: XTE J1550-564 feubota. k Makishimall2004j ) and 
4U 1630-47 (|Abe et al.l l2005ft. The spectral shape in the 
ASS is better represented by a slim disk model than the 
MCD model. If the spectra in this state are fitted with 
the MCD model, the following three anomali es appear 
(jKubota. k Makishimal l200i I Abe et~aTll2005l ): (i) Ex- 
cess emission is found below 5 keV and above 10 keV. 
(ii) The innermost disk temperature is higher than that 
in the HSS. (iii) The innermost disk radius tends to 
decrease with the increasing the disk temperature. In 
the HSS, in contrast , the innermost disk temperature is 
rather constant fe.g. lEbisawa et aLl ll994'). Some of these 
anomalies can be solved either by using the slim-disk 
model instead of the M CD model (|Kubota. k Makishimal 
120041: 1 Abe et all 120051 ). or by introducing spectral hard- 
ening to the MCD model as a function of accretion rate 
(|Davis et al.ll2006l : iMcClintock et al.ll2007l h Although it 
is yet unclear if the ASS should be considered to be an 
independent state or an extension of the HSS, we treat 
the four states (LHS, HSS, VHS, and ASS) as different 
states in this paper. 

Similarly to the Galactic BHBs, some ULXs are known 
to show state transitions in timescales of months to 
years between two differe nt spectral states. In the case 
of IC342 sources 1 and 2, IKubota et ail (|2001f ) dubbed 
them the "hard" and the "soft" state, which respectively 
has a spectrum with a PL and a convex shape. The latter 
is brighter than the former by from a factor of a few to 
one order of magnitude. However, the relation between 
the spectral states of ULXs and those of Galactic BHBs 
is still unclear. Observational clues will be obtained by 
comparing the state transitions, including their frequen- 
cies and time scales of transitions. 

ULXs are less well frequently monitored than Galac- 
tic BHBs, and their time variation is less understood. 
Obviously, increasing the samples of ULXs is a key to 
conduct statistical studies of state transitions. Nearby 
(<5 Mpc) galaxies and interacting galaxy systems at 
a moderate distance (5-10 Mpc) are two suitable lab- 
oratories for this purpose. ULXs in nearby galax- 
ies are bright in flux and can be examined in detail. 
About 15 ULXs were studi ed within 5 Mpc, w hich in- 
cludes NGC 1313 X-l and 2 dMizuno et alj |2007h. IC 342 



so urces 1 and 2 (IKubota et al.l 12001 ), Holmbergll X - 
1 (|Goad et all 120061). M81 X-9 (|Tsunoda et alj |2006), 
M82 X-l (iMiyawaki et all 120091). and NGC 5204 X-l 
(jRoberts et alj|2006f ) . 

We focus on the other laboratory; ULXs in interacting 
galaxy systems. They are more distant and fainter in 
flux than those in nearby galaxies. However, the average 
number of ULX per galaxy is much larger, enabling us to 
monitor a large number of samples at the same distance 
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simultaneously. The num ber of ULXs in nine inter acting 
systems at 20-100 Mpc (|Brassington et al.l 120071 ) aver- 
ages to be >6 ULXs system -1 . 

Within 10 Mpc, the M51 system and the NGC 4490 
and NGC 4485 system (hereafter called NGC 4490/85) 
host the l argest number of UL Xs. The former at 
~8.4 Mpc (jFeldmeier et al.lll997l) hosts nine, while the 
latter at ^8 Mpc (|de Vaucouleurs et al.1 119761 ) hosts 
eight. Both systems have been observed several times 
in X-rays and thus are suitable to study long-term vari- 
ations of multiple ULXs at a time. The X-ray data from 
these two systems have not been uniformly analyzed yet. 
In this paper, we apply physical models to the spectra of 
these ULXs and give interpretations to the observed long- 
term spectral changes, in comparison with the Galactic 
BHB state transitions. 

2. OBJECTS 
2.1. M51 

M51 is one of the closest interacting galaxy systems, 
which is composed of two spiral galaxies, NGC 5194 
and NGC 5195. The sizes of NGC 5194 and NGC 5195 
are 27 a nd 14 kpc from the isopho tal diameters 1U2 
and 5.'8 (|de Vaucouleurs et all 17991 ;. respectively. The 
total mass of the M 51 system is ~1.5xlO u M© 
(|Kuno. k Nakail fl997l) . and the mass of NGC 5194 is 
larger than tha t of NGC 5194 by a factor of a few 
(|Schweizerlll977t) . 

Radio observations ide ntified many H II regions 
(|van der Hulst et al.l Il988| ) , indicating that the system 
has an on-going star forming activity at a star formation 
rate of ~4 M Q yr" 1 (|Scoville et al.ll2001l) . Two (and pos- 
sibly another) core-collapse supernovae (SNe) recently 
occurred in the galaxy . 

M 51 h as been observed man y times in X-rays by the 
Einstein dPalumbo et allll985l ). ROSAT (Marston ^et^IJ 
19951 lEhle et all I1995D. ASCA dTerashima et al.| 119981 : 



Ptak et al.l 11999ft. BevvoSAX ([Fukazawa et al.l 120011) . 
XMM-Newton (IDewangan et al.| 120051). and Chandra 



(|Liu et al.ll200a iTerashima. fe WilsonN2004l ) observato- 
ries. The nu mber of X-ray sources increased from three 
by Einstein (jPalumbo et al.l fl985l ) to 113 by Chandra 

4 See |http: / / web.oapd.inaf.it/supern/cat / 1 for detail. 
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Fig. 1. — Optical and X-ray images of M51. (a) SDSS t hree color image. Red, g reen, and blue respectively represents the r-, g-, and 
M-band intensity. Pluses show positions of the nine ULXs in Dcwangan ct al. (2005); the two brightest sources, focused in this study, are 
shown in red and the remainders are in green, (b) Smoothed XMM-Newton MOS-1 image in the 0.5—8.0 keV band. Solid blue circles 
indicate the source extraction region. 
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Fig. 2.— Same with Figured] but for NGC 4490/85. Pluses show positions of the eight ULXs in lFridriksson et all j2008ft . The dashed 
circle in (b) indicates the common background region. 



(Terashima, & Wilson 2004), which includes nine ULXs 
(IDewangan et al1l2005l). We follow the nomenclature by 
iTerashima. fc Wilson! <|2004f > . Figure [1] (a) shows posi- 
tions of the ULXs on the optical image by the Sloan Dig- 
ital Sky Survey (SDSS) data. All ULX s are d istributed 
along galactic arms. ITerashima et al.l ()2006f ) reported 
using the Hubble Space Telescope data that four ULXs 
(sources-9, 37, 69, and 82) are located near the center or 
at the rim of star clusters. They also identified one or 



more optical counterparts to six ULXs. 

Previously publish e d X-ray st u dies (e.g. 
IDewangan et al.l l2005t ITerashima et al.l l2006h were 
based on four X-ray data sets among seven. We reduce 
all the available data including three unpublished sets. 

2.2. NGCU90/85 

NGC 4490 and NGC 4485, the former being a 
spiral and the latter being an irregular galaxy, 
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are interacting with each other. NGC 4490/85 is 
one of the closest interacting galaxies as well as 
M 5 1, belonging to the CV nil group of galax- 
ies (Ide Vaucouleurs et al.l 119761: IViallefond et alJ Il980t 
IClemens et all I1998L 119991 : IFridriksson et all I2008T I. 
The size and the ma s s of NGC 4490 are ~15 kpc 
([Clemens, fe Alexander! |200l) and -1.6xl0 10 M© 
( Viallefond et al.l I1980D . respectively, while those of 



NGC 4485 are ~ 5.6 kpc (IClemens. fc Alexanderl I2002T) 
and ~2xl0 9 M (jClemens et al.lll999f ). respectively. 

Radio observations identified many H II regions in 
NGC 4490, indicating an on-going star f orming activity 
at a star formation rate of ~4 . 7 Mp, yr _1 (Cle mens et al.l 
119991: IClemens. fe Alexanderl I20Q2H . Two core-collapse 
SNe recently occurred in the galaxy 4 . 

NGC 4490 / 8 5 was obs e rved several times with 
the ROSAT (|Read et al.l 119971: iRoberts. fc Warwick! 
120001). XMM-Newton, and Chandra o bservatories 
(jRoberts et aI1l2002t IFridriksson et all [2008). The num- 
ber of k nown X-ray sources incre ased from five by 
RO SAT (IRoberts. fc Warwick! I2000D to 38 by Chan- 
dra (|Fridriksson et al.l I2008T ) ■ which include eight ULXs 
(ULX-1 to ULX - 8). W e follow the nomenclature by 
Fridriksso n et al.l (|2008l ). Figure [5] (a) shows positions 
of the ULXs on the SDSS image. ULX-2 belongs to 
NGC 4485, while the others to NGC 4490. ULX-6 is close 
to but offset from t he dynamical center of the galaxy 
(|Roberts et alj|200l . and the other ULXs in NGC 4490 
are distributed along galactic arms. Six of them were 
observed by the Spitzer Space Telescope. Five sources 
(ULX-2, 3, 4, 6, and 8) are likely to be an accreting 
X-ray source based on the detection of some charac- 
teristic features of highly ionized species in the mid- 
infrared spectra, while the remai ning one (ULX-1) i s 
more likely to be an SN remnant (I Vazquez et al.ll2007l) . 
Nevertheless, X-ray flux variation was found from ULX-1 
(Fridriksso n et ai1l2008l ). leading to the speculation that 
this source is a ULX associated with the SN rem n ant. 

Previous X-r a y st udies (|Roberts et al.l 120021 : 
IFridriksson et al.1 l2008| ) presented long-term varia- 
tions in flux and color in several ULXs and gave some 
phenomenological analysis. We present long-term 
spectral variations of all the bright ULXs along with 
more physical models. 

3. OBSERVATIONS AND DATA REDUCTION 

We selected data sets with an exposure time longer 
than 10 ks for detailed spectral analysis. Throughout 
this paper, we use X-ray events in the 0.5-8.0 keV en- 
ergy band. For M51, we analyzed three Chandra and 
four XMM-Newton archived data sets (Table [J). We la- 
beled the observations as Cl-3 and Xl-4, resp e ctively . 



CI, C2, C3, and XI were u sed bv iLiu et all (2002) 



Dcwangan et al. 



2005) 



ITerashima. fc Wilson! (|200~ 
and iTerashima et al.1 (|2006[ ) . 

For NGC 4490/85, we reduced all the archived data 
sets of Chandra and XMM-Ne wton (Table ffl). The sets 
are the same with those used in [F ridrikss on et al.1 (|2008f ) 
and lGladstone. fc Roberts! (|2009h . We labeled the Chan- 
dra and the XMM-Newton observations as Cl'-3' and 
XI', respectively. 

3.1. Chandra 



Th e Chandra X-ray Observatory (|Wcisskop f et all 
12002ft has an unprecedented spatial resolution of ~0'.'5 
at the optical axis. The Chandra observations were 
conducted using the Advanced CCD for Imaging and 
Spectroscopy (ACIS; iGarmire et al.l 120031 ) . which cov- 
ers an energy range of 0.5-8.0 keV. The data were 
taken with the full frame mode with the aim point in 
the S3 chip. The field of view (FOV) of the chip is 
~500"x500", which contains the entire system of both 
M51 and NGC 4490/85. 

Using the Chandra Interactive Analysis of Observa- 
tions (CIAO) vers ion 4.0 and the ACIS Extract package 
(|Broos et alJl2002ft version 2008-06-06, we extracted the 
source and background events and constructed the en- 
ergy spectra. The source events were accumulated from 
a region around each source encircling 90% of photons 
of a point-like source. The background events were from 
an annulus around each source. 

We list the number of counts fo r the ULXs identi- 
fied in M51 (iDewangan et al.l 12005) and NGC 4490/85 
(|Fridriksson et al.l 12008) in Tables [5] and G2 respectively. 
Among all the ULXs, we focus on the sources with a to- 
tal count of more than 1000 at least in one observation, 
which are practically bright enough for our spectral anal- 
ysis. Seven sources are thus selected; sources-69 and 82 
for M 51, and ULX-2, 3, 4, 6, and 8 for NGC 4490/85. 

3.2. XMM-Newton 

The XMM-Newton (jJansen et al.ll2001l ) observatory is 
equipped with the European Photon Imaging Camera 
(EPIC), which is comprised of three X-ray CCD cameras. 
Two of the cameras are MOS arrays (MOS-1 and MOS-2) 
having an on- axis spatial reso lution of ~5" and s ensitiv- 
ity in the 0.15-12 keV band (jTurner et al.ll2001[) . while 
the remaining one is a pn array having an on- axis resolu- 
tion o f ~6" and sensitivity at 0.15-15 keV (jStriider et al.l 
2001). Both instruments have a FOV with a radius of 
30-40', and cover the two studied systems entirely (Fig- 
uresQJ) andEt"). In the M51 observations, the thin (XI, 
X3, and X4) and medium (X2) filters were used. All 
cameras were operated with the full frame mode. In the 
NGC 4490/85 observation (XI'), the medium filter was 
used. The two MOS and pn cameras were operated with 
the full frame mode and the extended full frame mode, 
respectively. 

We used the Science Analysis System (SAS) version 
7.1.0 for extracting events and generating response files. 
Some high background time intervals were seen in all ob- 
servations. For three observations (M51 X2, X3, and 
X4), we excluded intervals with count rates larger than 
the average by more than 3er. For the remaining obser- 
vations, we used all exposures because the effect is negli- 
gible with only a few counts in extracted source counts. 

We focus on the seven sources selected in the Chandra 
data (Tables [5] and [3]). The source extraction regions are 
shown in Figures [1] (b) and [2] (b). We chose the source 
region with a radius of 9-20" to avoid contamination 
from other sources. In M 51, the background regions were 
extracted from an annulus around each source to offset 
the global diffuse emission. In contrast, in NGC 4490/85, 
a common background region devoid of bright sources 
(Figure [5p) was adopted. The background events have 
~100 counts for MOS and ^300 counts for pn. 
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TABLE 2 

Source and Background Counts for M51. 
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ocatcd in a gap or 


a dead area of a 


chip. 











TABLE 3 

Source and Background Counts for NGC 4490/85. 



Name a 






Counts' 5 








Chandra 


XMM-Newton 






CT 


C2' 


C3' XI' 
MOS-1 MOS-2 


pn 



ULX-1 


239 (0.2) 


424 (1.1) 


305 (0.7) 








ULX-2 


1243 (0.8) 


1772 (9.9) 


1436 (0.9) 


609 (23.5) 


499 (21.0) 


403 (50.2) 


ULX-3 


566 (0.3) 


1148 (1.2) 


1209 (1.0) 


339 ( 4.8) 


309 ( 4.2) 


554 (10.2) 


ULX-4 


455 (0.4) 


1291 (1.9) 


1375 (1.3) 


365 ( 4.8) 


419 ( 4.2) 


775 (10.2) 


ULX-5 


30 (0.8) 


374 (1.7) 


561 (1.5) 








ULX-6 


473 (0.8) 


299 (1.9) 


1347 (2.3) 


444 ( 8.5) 


481 ( 7.5) 


897 (18.1) 


ULX-7 


(0.3) 


1 (0.7) 


749 (0.9) 








ULX-8 


780 (0.4) 


1333 (1.2) 


1762 (1.4) 


525 (23.5) 


484 (21.0) 


811 (50.2) 


a The nomc 


nclatures follow 


' liMdriksson et 


al.| 12003). The 


five brightest ULX 


s, focused in 


this paper, are 



shown in the bold font. 
b Same with Table[2] 

4. DATA ANALYSIS AND RESULTS 

Hereafter, we call each exposure of a source as a "sam- 
ple" . The total number of the samples is 34, since there 
are two sources with seven exposures for M51 and five 
sources with four exposures for NGC 4490/85. 

4.1. Light Curves 

Using the X-ray timing analysis package XRONOS ver- 
sion 5.2.1, we created the background-subtracted light 
curves of the seven ULXs. We binned the curves with 
1000 s bin -1 and fitted them with a constant flux model. 
For the XMM-Newton observations, we considered that 
the flux is variable when neither of the MOS and pn light 
curves are fitted with a constant model. 

As a result, we found significant variations (95% sig- 
nificance) only from one source (M51 source-69). This 
source showed a significant short-term variation in most 
samples (CI, C3, X2, X3, and X4). We confirmed the 



pe riodicity o f ~7000 s in the flux in CI as was reported 
bv lLiu et all ((20021) . 

4.2. Spectra 

Figure [3] shows the background-subtracted spectra 
of all the samples. All were binned with at least 
20 counts bin -1 . All spectra are featureless. Some 
sources show flux variation among different observations. 

4.2.1. Fitting Models 

We used the X-ray spectral fitting package XSPEC 
version 11.3.2 for the spectral analysis. For the 
interstellar extinction, we u sed the wabs model 
(M orrison. fc McCammonl fl983|) : for the hydrogen col- 
umn densities, we assumed the fixed Galactic extinc- 
tion of 1.57xl 20 cm" 2 for M51 or 1.78 xl0 20 cm- 2 for 
NGC 4490/85 (iDickev. fc Lockmanl [19901) and a thawed 
additional extinction for each source. For the XMM- 
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Fig. 3. — Time-averaged spectra of the seven ULXs. The top panels show the Chandra ACIS-S3 spectra by the three observations; CI, 
CI' (black), C2, C2' (red), and C3, C3' (green). The bottom panels show the XMM-Newton merged MOS spectra by the four (M51) or 
one (NGC 4490/85) observations; XI, XI' (black), X2 (red), X3 (green), and X4 (blue). 

Newton EPIC data, we fitted the merged MOS and pn 
spectra simultaneously. 

We first applied a PL model and an MCD model 
(diskbb in XSPEC; lM~itsuda et al.lll984f ). which are com- 
monly used to fit the continuum emission of Galactic 
BHBs and ULXs. The MCD model has the innermost 
disk temperature T- ln as one of the parameters. Using the 
relation L^ Q \ = iir(Ri D /£,) 2 a(Ti n/ k) 4 , we can calculat e 
the innermost disk radius i?i n (Makishi ma et all [2000). 
Here Lboi is the bolometric luminosity, a is the Stcfan- 
Boltz mann constant, k ~ 1.7 is th e spectral hardening 
factor (Shimura, & Takahara 1995), which is the ratio of 
the color temperature to the effective temperature, and 
£ = 0.412 is t he correction factor for the inner bound- 
ary condition (|Kubota et al.l 1 1998T) . For non-spinning 
BHs, the BH mass -Mbh can be estimat ed from Rj n using 
the r elation i?; n = 8.86(A/bh/Mq) km (jMakishima et all 
I2000D . For the MCD model with kT in = 0.7-2.0 keV, the 
0.5-8.0 keV band contains >80% of the entire spectrum, 
thus we approximate the X-ray luminosity in this band 
as the bolometric luminosity. 

We also fitted all t he spectra with a slim disk model 
bv iKawaguchil (j2003l ) 5 . The "Kawaguchi model" takes 



into account of the Comptonization and the relativistic 
effects, and can constrain the BH mass and mass ac- 
cretion rate as model parameters. We fixed the viscous 

5 See http://heasarc.nasa.gov/xanadu/xspec/models/slimdisk.html 
for detail! 



parameter a to be 0.1 ()Vierdavanti et all l2006f) . This 
model was success fully applied to some bright ULXs, in- 
cluding M33 X-8 (|Foschini et a l.l l200l. NGC 1313 X-2 
NGC 4559 X-7, X-10 NGC 5204 X-l (IVierdavanti etahl 
12051 . and M82 X-l (lOkaiima et al.ll2006D . 

4.2.2. Fitting Results 

We summarize all the best-fit parameters in Ta- 
bles 1 and [5] for samples in M51 and NGC 4490/85, 
respectively. N^, r, and dM/dt indicate the absorp- 
tion column density additional to that in our Galaxy, 
the photon index, and the mass accretion rate, re- 
spectively. Using the absorption-corrected flux (,/xpli 
/x,mod> and /x,Kaw respectively for the PL, MCD, 
and Kawaguchi models), the luminosity is calculated as 
Lx,pl = 47tD 2 /^ pl , L x ,mcd = 27rD 2 / XMCD (cosj)~ 1 , 
and £x,Kaw = 27rZ? 2 / x Kaw (cosi) _1 . Here, D is the dis- 
tance to each galaxy, and i is the inclination of the disk. 
We assumed isotropic emission for the PL, and a mod- 
erate inclination for the MCD and Kawaguchi models 
(i = 45°). Errors of all the parameters are calculated at 
the 90% confidence level. Figure 2] shows examples of the 
spectra and the best-fit models. 

Although short-term flux and spectral hardness varia- 
tion are found in M 51 source-69 in CI, we averaged each 
exposure because the paucity of the counts do not al- 
low us to conduct time-sliced spectroscopy. This source 
was very faint in C2, and is fixed at 6.0xl0 20 cm" 1 , 
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TABLE 4 
Best-fit Parameters for M51. 



Name Data 
label 


Model* 


(10 22 cm- 2 ) 


r 


kT in 
(keV) 


^in 

(km) 


(Mq) 


dM/dt 
(L Edd c~ 2 ) 


(10 38 ergs s" 1 ) 


Red-x 2 
(d.o.f.) 


Source-69 CI 


P 

M 
K 


< 0.08 

< 0.03 

T>+ - 07 
u - io -0.07 


J - 1 -0.1 


9 q+1.2 
z -°— 0.4 


26± 2 | 




looj 


513 
71 


27 ±3 
18+ 3 
16+. I 


1.12(18) 
1.33(18) 
1.99(18) 


C2 


P 

M 
K 


0.06(fixed) 
0.06(fixed) 
0.06(fixed) 


3 1+ 4 ' 1 




19 ,+143 


io± 3 6 


< 2 


< 0.8 

< 0.4 

< 0.6 


0.63(1) 
0.14(1) 
0.18(1) 


C3 


P 

M 
K 


05 +u - u4 

u - ua -0.04 

< 0.01 
07+ 03 

u - u '-0.03 


1 5 +ul 

L - a -o.i 




60+ 16 
ou -10 


7+ 4 

'- 3 


iool 3 i 


29 ±2 
16 ± 1 
19 ± 1 


0.80(34) 
1.94(34) 
1.13(34) 


XI 


P 

M 
K 


14+ u - ii 

< 0.03 
u - uo -0.05 


L - y -0.2 


i 9+O.2 


ec+21 
°°-13 


< 9 


> 425 


12 ± 1 
6.1 ± 0.6 
7.5 ±0.7 


0.86(40) 
1.19(40) 
0.91(40) 


X2 


P 
M 
K 


n nq +u - u,:! 

U.UM_q 03 

0.00 


-1 7 +u.u 
i -'-o.i 


1 5+ ' 1 


sit 1 ! 


\2+_ 2 


lOOl 


254 
41 


52 ± 2 
30 ±1 
34 ±1 


1.18(110) 
2.03(110) 
1.33(110) 


X3 


P 

M 
K 


o 09 +u - UJ 

u ' uy -0.03 

0.00 

10+ 002 
u ' iu -0.02 


1 fi +ul 
^•"-o.i 


1 7 +0.i 
L -' -0.1 


581 8 9 


10t \ 


100+ 


154 
34 


41 ±2 
25 ± 1 
27+. 2 


1.58(113) 
2.23(113) 
1.76(113) 


X4 


P 

M 
K 


u ' iu -0.03 

< 0.01 

13+ 002 
u - lo -0.02 


1 5+"- 1 


1 S + ' 2 


60t 10 7 


13± 2 


ioot 23 


58 ±2 
35 ± 2 
381 I 


1.41(111) 
1.80(111) 
1.71(111) 


Source-82 CI 


P 

M 
K 


o n+ u - Ub 

u - L1 -0.06 

< 0.01 
01+ 05 


9 o + U.3 
z -°-0.2 


7+ ai 
u -'-o.i 


222+Jji 


34t 7 


7+ 


3 
1 


28 ±3 
13 ± 1 
15± I 


0.84(30) 
1.48(30) 
0.95(30) 


C2 


P 

M 
K 


02 +u - ud 

u - uz -0.02 

< 0.01 

< 0.01 


1 Q+U 


9+ 01 
u ' y -0.1 


-.91+38 

izl -36 


121 3 




23 
5 


16 ± 1 
8.5 ±0.9 
11 ± 1 


1.13(30) 
2.56(30) 
1.31(30) 


C3 


P 

M 
K 


u ' iz -0.04 

< 0.01 
04+ 03 


z ' i -0.1 


9+ 01 
u - J -0.1 


1 ic+21 
146 -20 


221 1 " 


9i 


7 
2 


25 ± 2 
12 ± 1 
14 ±1 


1.60(34) 
2.27(34) 
1.72(34) 


XI 


P 
M 
K 


21 +u - u ' 
< 0.03 

.08±r i 


2 6+ u - 3 

z -°-0.2 


7+ ' 1 
u - '-0.1 


lq1 +39 


37+. ? 


si 


1 
1 


271 S 
11 ± 1 

13 ± 1 


1.18(47) 
1.91(47) 
1.45(47) 


X2 


P 

M 
K 


< 0.04 

< 0.01 

< 0.01 






107+J 4 


11+ g 


20j 


58 
8 


19 ±2 
11 ± 1 
13 ± 1 


0.97(31) 
2.58(31) 
1.22(31) 


X3 


P 

M 
K 


o n+ u - u4 

< 0.01 

n m+ 03 

u - ui -0.01 


n n+U.2 

z - z -o.i 


o-s±8:i 


152±H 


171 5 3 


7+ 


4 
2 


17 ±1 
8.1 ± 0.5 
9.6 ±0.6 


0.87(75) 
1.80(75) 
1.13(75) 


X4 


P 

M 
K 


14+ U - Ub 

< 0.01 
04+ 04 


n 9+U.I 

z - z -0.2 


0.9+83 


132±|g 


19l ? 


si 


7 
2 


20 ±2 
9.6 ±0.6 
11 ± 1 


0.83(62) 
1.47(62) 
1.03(62) 



The abbreviations for the models: "P" for PL, "M" for MCD (the diskbb model in XSPEC), and "K" for Kawaguchi model. 



which is the median value obtained from the other six 
observations. 

5. DISCUSSION 

5.1. ULX Spectral States and Transitions 

Based on the spectral fitting results presented in § 4, 
we investigate the spectral states of ULXs and their tran- 
sitions. It is desirable that the states be defined for in- 
dividual sources. However, given the low photon counts 
of each source, starting with individual details may lead 
to an unclear view. We thus first investigate the over- 
all trend of all samples to constrain states of ULXs as a 
whole. In this section, we develop a discussion under the 
working hypothesis that all ULXs have a similar mass 



and have the same states and follow a similar pattern 
of state transitions. The working hypothesis is checked 
with individual sources in § 5.2. 

We constructed a histogram of the samples fitted suc- 
cessfully with the PL or the MCD model that represents 
PL-like and curved spectra, respectively (Figure [5]). We 
consider that the fitting with a null hypothesis probabil- 
ity of >5% is successful. The best-fit luminosities derived 
by the PL fitting are used, because the derived luminosi- 
ties hardly depend on the fitted models. 

In the plot, it is noticeable that the fraction of suc- 
cessful fits by the MCD model is higher within the range 
of 3-6 xlO 39 ergs s _1 for 20 samples and that by the PL 
model is higher within the range of 1.5-3xl0 39 ergs s _1 
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TABLE 5 

Best-fit Parameters for NGC 4490/85. 



Name 


Data 
label 


Model a 


(10 22 cm" 2 ) 


r 


kT- ln 
(keV) 


^in 

(km) 


M BH 
(Mq) 


dM/dt 
(^Edd c~ 2 ) 


(10 38 ergs s _1 ) 


Rcd-x 2 
(d.o.f.) 


ULX-2 


CI' 


P 
M 
K 


n 21 + aos 

u - zi -0.07 

o n c i +0 ' 05 

u ' uo -0.05 
n 99 +0.05 
u - zz -0.04 


1 -°-o.a 


1 6+°- 2 


7 q+22 


12l 12 


ioolf° 6 


52 ±4 
29 ± 2 
35l 3 


1.69(28) 
1.34(28) 
1.68(28) 




C2' 


P 
M 
K 


o 21+ U - Ub 

u - zl -0.06 

u ' uo -0.04 
n -,0+0.04 


1 -°-0.2 


i o+0.2 
A - a -0.1 


Q9+22 
az -18 


18lJ° 


> 13 


40 ±2 
20 ±1 

251 3 


0.94(35) 
0.77(35) 
0.87(35) 




C3' 


P 
M 
K 


0.28+S 

o.o8ig : « 

n 24+ ' 06 


1 R+0' 1 
L -°-0.1 


0.1 


79+18 
,z -15 




> 14 


34 ±2 
18 ± 1 
22 ± 2 


1.63(28) 
1.28(28) 
1.50(28) 




XI' 


P 
M 
K 


o 2 7+ u ' u/ 

u - z '-0.07 
u ' uz -0.02 

n ifi+ ' 04 

u - iD -0.04 


r, 1+U.l 

z - i -0.2 


1 2+ - 1 


97+21 
9 '-20 


231° 


lei 1 " 


401 t 
19 ± 1 
24 ±2 


1.36(48) 
1.60(48) 
1.40(48) 


ULX-3 


CI' 


P 
M 
K 


S7+ 24 

u - il -0.21 
n 65+ - 21 


2 0+°' 4 
z - u -0.3 


1 5+°- 4 

- 1 - — 0.3 


65lg 


17+16 
1 ' -13 


> 12 


36l i2 
18 ±2 
22 ±4 


1.05(22) 
1.14(22) 
1.06(22) 




C2' 


P 
M 
K 


1 QQ+U-^b 

l.OM_ 23 

0.791°^ 


2 4+0.3 
Z ' 4 -0.3 


i 9+O.2 

J -- z -o.i 


113l 3 ° 


40l | 


14+ % 


581- 
27 ±2 
33l % 


0.63(30) 
0.67(30) 
0.60(30) 




C3' 


P 
M 
K 


Los!":!? 

0.601°;" 
0.891°;" 


9 -i+U.2 

z -0.2 


1 4+0-2 

-"■-^-o.i 


rst\ 9 6 


29ll0 


171 1 ! 


Bit 7 
23 ± 2 
29l I 


0.90(33) 
0.84(33) 
0.82(33) 




XI' 


P 
M 
K 


U.80_ Q 20 

39+ 014 
u - ' 4 -0.10 


1 c+U-2 
1 -°-0.1 


1 s+0.3 
1 - 8 -0.2 


501- 


121" 


> 23 


441 I 
25 ± 2 
291* 


1.00(33) 
1.10(33) 
1.00(33) 


ULX-4 


CI' 


P 
M 
K 


49+P u 

22+ 13 
u ' zz -0.12 

41 +°' 13 
u '^ 1 -0.12 


J - 8 -0.2 


1 5+ ' 4 
J -°-0.3 


55+ 29 


12l 6 9 


> 12 


25+ i 
13+ I 
161 1 


1.08(17) 
1.15(17) 
1.08(17) 




C2' 


P 
M 
K 


°- 63 -0.15 

29+ 11 
u ' zy -0.10 

54+° 11 


1 q+O-i 


1 5+ ' 2 
0.2 


81+« 


26+ 9 
ZD -20 


> 16 


27 ±2 

33l | 


1.36(29) 
1.24(29) 
1.29(29) 




C3' 


P 
M 
K 


64+°' 13 
u -°+-0.06 

u " 54 -0.08 

0.571°;?? 


1 g+0. -i 


i 5+°- 2 


70111 


wis 


> 16 


42 ±4 
22 ± 1 
27 ±2 


1.63(28) 
1.42(28) 
1.53(28) 




XI' 


P 
M 
K 


0.601°;- 
o 3n+° 08 

u - au -0.07 


2 0+ a2 
z - u -0.1 


1 Q + 0.1 


99l?? 


341 5 8 


1711 


501^ 
26 ±2 
311 3 


1.32(44) 
0.98(44) 
1.14(44) 


ULX-6 


CI' 


P 
Al 
K 


3-5+0' 23 

0,l0l°;« 
0-28±g;« 


1 S + 0-^ 
J--8-0.3 


1 4+0- 5 
- L - 4 -0.3 


°°-22 


io+ ; 


> 11 


23+ « 
12l 2 
15 ±2 


0.72(17) 
0.79(17) 
0.73(17) 




C2' 


P 
M 
K 


u.ou_ 20 
0.35l°; 2 ? 


2 6+ ab 
z - D -0.4 


u.y_ 2 


64l« 


12l 6 


is! i 


7+ 7 
'- 3 
o n +0.6 

o s+1.3 
3-°-1.0 


1.09(10) 
1.45(10) 
1.24(10) 




C3' 


P 
M 
K 


0.88l°;f 4 

45+° 10 
u ' 4o -0.09 

65+ 010 


2 4+°' 2 
z -^-0.2 


1 1+ - 1 


1241H 


37l | 


ill I 


53 +1 ° 
2ll 2 
26 ±2 


1.49(31) 
0.93(31) 
1.06(31) 




XI' 


P 
M 
K 


u - 41 -0.08 

11 + 005 
U - il -0.05 

U.OU_ 05 


o i+U.l 
z ' i -0.1 


1 4+0-1 
1 - <t — 0.1 


86l« 


27l 7 


171 I 


44 ±4 
23 ± 1 
27l % 


1.12(52) 
1.07(52) 
1.02(52) 


ULX-8 


CI' 


P 
Al 
K 


U - 8b -0.20 

0.461°;^ 
0.68l°« 


9 o+U-3 
z -°-0.3 


i 9+O.2 


H6±l! 


37l 8 


12l | 


56l ib 9 
24 ±2 
30 ±4 


1.45(20) 
0.96(20) 
1.13(20) 




C2' 


P 
Al 
K 


0.901°- 

43+ 010 
u ' 4J -0.09 

60+ 12 
u - ou -0.09 


2 6+ a2 
z - D -0.2 


1 0+ - 1 

1 - u -0.1 


14711 8 


40l 13 


9+ i 

a - 2 


46li 
19 ±2 
23+ 3 


1.28(26) 
0.75(26) 
0.82(26) 




C3' 


P 
Al 
K 


0.74l°- 

0.38l°;°? 

Q61 +0.08 
u ' D1 -0.07 


2 0+ a2 
ZiU -0.1 


1 4+0-1 

-■-*-o.i 


90+1S 


32+ 4 
oz -10 


20l 12 


59l I 
28 ± 2 
35 ±4 


1.52(36) 
1.47(36) 
1.41(36) 




XI' 


P 
Al 
K 


064 +u.io 

u -°*-0.09 

26+ a07 
u ' ZD -0.06 
sg+0-08 


9 r+U.2 
z '°-0.2 


1 0+ - 1 


152l 2 2 8 5 


381 10 , 


iol J 


54l 9 
21 ± 1 
26l 2 


0.88(58) 
0.92(58) 
0.86(58) 



a Same with Tabic g] 
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Fig. 4. — Examples of spectral fitting in M51 (a— d) and NGC 4490/85 (e— h). Spectra are shown for different sources with different 
instruments as indicated by the labels and colors. The top panels show the spectra (pluses) and the best-fit PL model (solid histogram), 
and the other panels show the residuals from the best-fit by a PL, MCD, and Kawaguchi model. 
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10 20 50 100 

Lx.pl (10 38 ergs s" 1 ) 

Fig. 5. — Luminosity distribution of (a) the number of samples 
and (b) the fraction of successful fits by the PL (black) and the 
MCD (red) models. The 1<t uncertainty, which is calculated from 
the Poisson distribution, is given in (a). 

for 11 samples (Figure \5jp) . In order to test the claim 
statistically, we simulated numerous spectra containing 
the two components (PL and MCD) with randomly cho- 
sen parameters with varying contrasts, fitted them with 
one component model of either the PL or the MCD, and 
derived the fraction of the successful fits. As a result, we 
found no trend that the MCD model is favored toward 
brighter luminosity. Indeed, the PL model was favored. 
Thus, we consider that the higher fraction of successful 
MCD fits toward the brighter luminosity is not a statis- 
tical artifact. 

The bimodal structure in the luminosity is unlikely 
to be explained by the bimodal distribution of the BH 
mass or the inclination of the system. This is because 
the distribution of the mass or the angle should be con- 
tinuous, thus such a bimodal distribution is not natu- 
ral. Rather, we speculate that there are two distinc- 
tive states corresponding to each luminosity range. We 
hereafter call them the "bright' state and the "faint" 
state. The two states match the "so ft " and the "hard" 
states of the ULXs in iKubota et al.l ()2001[ ). Based on 
the number of samples in each state, the frequency ratio 
bright: faint is about 2:1. The transition between the two 
occurs at ~3xl0 39 ergs s -1 (hereafter the "border lumi- 
nosity" ) . iGladstone. fc Roberts! (|2009t ) also claimed that 
the changes in spectral state occur at a similar luminosity 
of ~2xl0 39 ergs s" 1 using samples in NGC 4490/85. 

Some samples do not belong to either of the two 
states. Three samples (M51 source-69 in C2, XI, and 
NGC 4490/85 ULX-6 in C2') show luminosities fainter 
than the tail of the fainter range. These spectra may 
belong to the third state, which we call the "dim" state. 
We cannot constrain its representative spectral model for 
the paucity of statistics. 

Two types of luminosity variation are recognized. We 
designate the variations within one state as the "intra- 
state" variations, and those across two or more states as 
the 11 inter- state" transitions. 

5.2. Spectral Variations of Individual Sources 

Now, we look into details of individual sources. The 
following result shows that the BH mass is in a range of 
30-40 Mq for sources which we can estimate a mass with 
a representative model. This supports our assumption 



to derive the overall trend in § 5.1. All the samples can 
be considered to be in one of the three states defined 
in § 5.1. The spectra in the bright state is explained 
better with a curved model and those in the faint state 
with a PL model. For the curved model, we see that 
the Kawaguchi model better explains the sources with 
sufficient constraints. 

For each individual source, we discuss the states of 
all samples and the most likely representative model for 
each state. The result of spectral fitting is employed, 
which is mainly used to distinguish PL-like or curved 
spectra. If both the MCD and the Kawaguchi models 
give an acceptable fit for a curved spectrum, we further 
use a plot of Lx.mcd against kT ln (Figure [5]) to distin- 
guish the two. The plot is commonly used for Galactic 
BHBs and ULXs. The two parameters are correlated in 
a different manner in two different physical conditions: 
standard disk or slim disk. If a system has a standard 
disk (the MCD model), the relation Lx,mcd °c T; n is 
obtained, which is interpreted that the innermost disk 
radius remain s constant against varying accretion rates 
in th e state (IShakura. fc Sunyaevi 119731 : lEbisawa et al.l 
119931: iGierlihski. fc Pond 120040 ^ a system has a slim 
disk (the Kawaguchi model), the relation Lx.mcd °c Tj n ' 3 
is obtained, where f3 < 4 ijWatarai et al.l l2000h . The 
powe r decreases from 4 as an increasing mass accretion 
rate (jWatarai et al.|[200lT) . For samples with a successful 
fit by the MCD or Kawaguchi model or both, we derived 
the best-fit Mbh value and examined if it is consistent 
among all samples belonging to a source. 

5.2.1. NGC4490/85 ULXS 

We consider that ULX-8 exhibited an intra-state varia- 
tion within the bright state. This is because the luminos- 
ity variation (4.6-6.0 x 10 39 ergs s" 1 ) happens only above 
the border luminosity. 

We speculate that this ULX is in the slim disk state 
based on the two lines of evidence. The first is that 
the scatter in the £x,MCD - Ti n plot (Figure [5^) follows 



^X,MCD 



cx TiJ, where = 1.0 



+0.7 
-0.4- 



The value is consis- 



tent with the slim disk (J3 < 4) and not with the standard 
disk (/3 ~ 4). We note that the value is not an apparent 
slope caused by an additional spectral component (e.g. 
an additional PL component generally found in the HSS), 
because we obtained an even flatter slope {/3 = 0.6t a 7 3 ) 
when the spectra were fitted by two component models 
(the PL plus the MCD). Here, we fixed the photon index 
of the additional PL component to be 2.0 for the lack of 
statistics to constrain the value. 

The second evidence is that a BH mass consistent 
among all the samples is found for the Kawaguchi 
model, but not for the MCD model. Figure [7] shows 
the best-fit parameters using the MCD and Kawaguchi 
models. In the MCD fitting (Figure [7^,), a value 
(i? in /km) 2 /(D/10 kpc) 2 (cosi)" 1 cx M BH 2 consistent 
among the four samples was not found within 90% sta- 
tistical uncertainty, which is unphysical. In other words, 
although spectra of ULXs in the bright s tate are com- 
monl y analyzed using the MCD model fe.g. lKubota et all 
2001), it now turns out that the MCD model does not 
describe the long-term spectral variations of the bright 
state. On the other hand, in the Kawaguchi model fit- 
ting (Figure [TJd) , the range Mbh = 35-38 M© was found 
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Fig. 6. — Plot of X-ray luminosity Lx.MCD against the innermost disk temperature kTi n based on the MCD fitting results. The results of 
different observations are shown in different colors. The dashed lines show the Lx,MCD °c Ti n 4 relation with several representative masses, 
while the dotted line indicates the Eddington luminosity for the standard disk. 
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Fig. 7. — Significance contours of the best-fit physical parameters by two different models for NGC 4490/85 ULX-8; (a) fcT; n versus the 
normalization (ocAfBH 2 ) by the MCD model, (b) dM/dt versus Mbh by the Kawaguchi model. The results in different observations are 
shown with different colors: Cf (black), C2' (red), C3' (green), and XI' (blue). The two confidence levels are shown for 68% (thick solid) 
and 90% (thin solid) ranges for the two parameters. The inset at the top right in (b) is a close-up view of the center in the linear scale. The 
dark and light gray shaded ranges indicate the common A/bh ranges among the four samples with a 68% and 90% confidence, respectively. 



to be consistent within 90% among the four samples. 
Thus, at least for this source, the Kawaguchi model is a 
more favorable description for its spectral variation. 

Now we found a common BH mass range acceptable for 
the four samples, we fitted the four spectra by tying the 
BH mass to obtain a more stringent constraint. The best- 
fit value of M B h is 37 ±2 Mq (Table [6]). Resultantly, the 
Eddington ratio (Lx/^Edd) is 0.41-0.63, which is similar 
to those of other ULX spectra considered to be in the slim 
disk state (0.36 for NGC 1313 X-2 and 0.52 for NGC 4559 
X-7; IVierdayanti et al.ll2006| ) . 

5.2.2. M51 Source-82 

Source-82 shows variability among samples, which we 
consider to be an intra-state origin within the faint 
state. This is because the luminosity variation (1.6- 
2.9 xlO 39 ergs s _1 ) happens only below the border lu- 



TABLE 6 

Best-fit parameters of the Kawaguchi model 
for ULX-8. 



Data dM/dt L x ,Kaw 

label (fO 22 cm~ 2 ) (L Edd r 2 ) (10 38 ergs s" 1 ) 



Cf 0.68I O o9 12+j 30 ±3 

C2' 0.55+° :° l 10tJ 22± I 

C3' 0.68jlg;S5e 15t? 36 ± 2 

XI' 0.38±g : g| 11±{ 26± | 

A/ BH a (Mq) : 37 ±2 

Red-x 2 (d.o.f.) 1.04(143) 



a The BH mass is fixed among the four independent fit- 
tings. 

minosity. 

We consider that the source stayed in the state rep- 
resented by the PL model in all samples based on the 
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Fig. 8. — Significance contours of the best-fit physical parameters 
(dM/dt versus Mbh) by the Kawaguchi model for M 51 source-82; 
The results in different observations are shown with different colors: 
CI (black), C2 (red), C3 (green), XI (blue), X2 (gray), X3 (light 
red), and X4 (light green). The two confidence levels are shown for 
68% (thick solid) and 90% (thin solid) ranges. 

following three reasons: (i) The majority of samples (six 
out of all seven) are successfully explained by a PL model, 
(ii) None of the samples is explained by an MCD model, 
(hi) Although the Kawaguchi model reproduces all the 
samples, the model is unlikely because we do not find 
Mbh consistent among all the samples (Figure [8]). 

In one sample (XI), source-82 shows a steep photon 
index of 2. (Table SJ. In contrast, in the other six 
samples, the derived index (r ~ 1.8-2.2; Table E) is flat- 
ter than that in XI. 

5.2.3. M51 Source-69 

Source-69 clearly caused an inter-state transition. This 
ULX exhibits the largest luminosity variation by at least 
a factor of 70, and shows a transition among the bright 
state (X2, X3, and X4; 4.0-5.8xl0 39 ergs s" 1 ), the faint 
state (CI and C3; 2.6-2.9 xlO 39 ergs s _1 ), and the dim 
state (C2 and XI; <1.2xl0 39 ergs s^ 1 ). 

We discuss the representative model for the three 
states of the ULX. First, in the bright state, none of the 
PL, MCD, and Kawaguchi models explains all the three 
samples (Table 0]). In particular, two samples (X3 and 
X4) were not reproduced by any models. We employed 
two component models (the PL plus the MCD or the 
Kawaguchi model), which did not improve the fit. The 
representative model was not found for this source in this 
state. 

Second, the PL is the representative model for the faint 
state, because it is the only model to reproduce spectra 
of all samples. The derived index is quite flat of r < 1.5. 

Finally, for the dim state, any models can reproduce 
the spectra and the representative model was not con- 
strained. We note, however, that the spectra in this 
state is softer (r ~ 1.9-3.1) than in the other two states 
(r < 1.7). 

5.2.4. NGCU90/85 ULX-6 

We consider that ULX-6 caused an inter-state tran- 
sition because the variation spans the three different 
luminosity ranges as M 51 source-69: the bright state 
(C3' and XI'; 4.4-5. 3xl0 39 ergs s" 1 ), the faint state 
(CI'; -2.3xl0 39 ergs s^ 1 ), and the dim state (C2'; 
-7.2xl0 38 ergs s" 1 ). 



For the bright state, both the MCD and the Kawaguchi 
models can explain all the spectra, and yielded a BH 
mass consistent among the two samples. However, the 
flat slope 08 = OAtfy of the L X , M CD-T in plot (Fig- 
ure^) favors the slim disk interpretation (f3 < 4). Thus, 
the representative model of the ULX in this state is ex- 
pected to be the Kawaguchi model. We tried a simultane- 
ous fitting by tying the Mbh value using the Kawaguchi 
model, and obtained the best-fit BH mass of 33±5 M Q . 
The Eddington ratio is 0.53-0.55. 

For the faint and dim states, all the models can explain 
all the samples, thus the representative model was not 
constrained. In the dim state, the PL fitting yielded the 
best-fit power of ~ 2.6, which is similar to that found in 
the dim state of M51 source-69. 

5.2.5. NGC 4490/85 the other ULXs 

The remaining three sources have no clear features, 
thus the representative model and the BH mass can 
hardly be constrained. We briefly describe these sources. 

Both ULX-2 and ULX-3 show intra-state varia- 
tion as their luminosity (Lx.pl ~ 3.3-5.2 and 3.5- 
5.8xl0 39 ergs s _1 , respectively) stays in the bright state. 
For ULX-2, no model explains all the samples, thus the 
representative model was not found for this source. For 
ULX-3, any models can reproduce all the spectra. The 
Lx,MCD~Ti n plot does not separate the slim disk or the 
standard disk state, because the data cannot be fitted 
with a linear relation. A BH mass consistent among 
all the samples is found both for the MCD and the 
Kawaguchi models as 8 ± 1 Mq and 28^y M Q , respec- 
tively. 

ULX-4, on the other hand, is considered to have ex- 
hibited an inter-state transition. The luminosities of the 
three samples are in the bright state (C2', C3', and XI'; 
4.2-5. lxlO 39 ergs s" 1 ), while that of the other is in the 
faint state (CI'; ~2.5xl0 39 ergs s^ 1 ). In both states, any 
models explain most of the samples, and the representa- 
tive model cannot be constrained for the same reasons 
as ULX-3. If we estimate the BH mass using a simulta- 
neous fitting of the two samples in the bright state, the 
BH mass is obtained as 9 ± 1 M (MCD) and 26± 1 q M Q 
(Kawaguchi). 

5.3. Relations with Galactic BHBs 



Many previous s t udies (e.g. iKubqta et all 



20011: iDewanean et al. 


120051: Winter et al. 


2006; 


Vierdavanti et all 12006 


; iGladstone, & Roberts! 


2009; 



Godet et al. 2009) suggested that the states of ULXs 



are related to those of Galactic BHBs. In the present 
study, we uniformly analyzed 34 spectra of seven ULXs 
in two interacting galaxies. As shown in § 5.2, we found 
some patterns of luminosity and spectral variations. We 
proceed to compare each state of the ULXs with those 
of Galactic BHBs. 

For the bright state, especially for NGC 4490/85 ULX- 
8 and 6, we speculate that this state corresponds to 
the ASS based on the following two reasons: (i) The 
two sources with the slim disk are the most likely in- 
terpretation (§ 5.2.1 and 5.2.4). (ii) The Eddington ra- 
tios (0.41-0.63 and 0.53-0.55) are similar to those of two 
Galactic B HBs considered to be in the A SS: XTE J1550- 
564 (~0.4: lKubota. fc Makishimal[200l and 4U 1630-47 
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fO.25-0.35: [Abe et all [20051) . However, the bright state 
of M51 source-69 and NGC 4490/85 ULX-2 do not show 
such features of the ASS. This fact suggests that the 
boundary luminosity to separate the bright and the faint 
states may be correct overall, but varies source to source. 

For the faint state, the representative model is PL for 
sources which we can determine a model (M 51 source-82 
and 69). The VHS and the LHS in Galactic BHBs have 
a PL shaped spectrum with the photon index r of ^2.5 
and 1.5-1.9, respectively. However, our samples in the 
faint state show a wide PL index range of r ~ 1.1-2.6. 
Because a state with such a steep index of 1.1 has not 
been known so far, the faint state may be different from 
the well-known state in Galactic BHBs. 

Finally, we briefly mention the dim state. We consider 
the dim state seen in some sources is distinctive from the 
faint states. This is because this state has lower luminos- 
ity and steeper spectral index, when fitted with PL, than 
the faint state. The photon statistics is too low to con- 
strain the representative model, and their counterpart 
state in the Galactic BHBs is unknown. 

6. SUMMARY AND CONCLUSIONS 

We have analyzed X-ray spectra (34 samples) of seven 
bright ULXs in the interacting galaxy systems M51 
and NGC 4490/85 using the archived data from multi- 
ple Chandra and XMM-Newton observations. We con- 
structed a histogram of luminosities, in which we found 
a hint of three distinctive states: (i) The bright state 
with the brighter luminosity range (3-6 xlO 39 ergs s _1 ) 
for 20 samples, (ii) the faint state with the fainter lumi- 
nosity range (1.5-3xl0 39 ergs s -1 ) for 11 samples, and 
(iii) the dim state with luminosities below the fainter 
range (<1.2xl0 39 ergs s _1 ) for three samples. On the 
whole, the spectra in the bright state are represented by 
the MCD model (standard disk) or the Kawaguchi model 
(slim disk), while those in the faint state are represented 
by the PL model. We classified the flux changes of all 
sources into the intra-state variation and the inter-state 
transition. Four ULXs (M51 source-82, NGC 4490/85 
ULX-2, 3, and 8) exhibited the former, while the remain- 
ing three (M51 source-69, NGC 4490/85 ULX-4, and 6) 
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exhibited the latter. 

All the spectra of two ULXs (NGC 4490/85 ULX-8 and 
6) in the bright state can be reproduced by a slim disk 
model by Kawaguchi (2003) with a constant BH mass 
of 37 ± 2 and 33^ M Q , respectively. In particular, for 
ULX-8, a common BH mass was not found with the MCD 
model fitting. This suggests that the slim disk model is 
a more appropriate explanation at least in some cases of 
the " bright" state, in which ULXs harbor BHs of a mass 
of <40 M Q . 

From the results of spectral fitting, we compare each 
state of the ULXs with those of Galactic BHBs. We pro- 
pose that the bright state of two ULXs (NGC 4490/85 
ULX-8 and 6) corresponds to the ASS of Galactic BHBs. 
For the faint state, the representative model is PL. How- 
ever, it is different in the range of photon index from the 
VHS of Galactic BHBs, which also shows a PL spectrum. 
We thus consider that the faint state is unique to ULXs. 
The nature of the dim state is unconstrained, but has 
the most steep spectrum. 
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